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COMPRESSION CONTEXTUELLE D'IMAGES NUMERIQUES 



L'invention concerne le domains de la compression/decompression de 
5 donnees d'images numeriques par transformee en ondelettes. 

Le codage par transform^ en ondelettes consiste classiquement a 
coder une image en trois §tapes. Une premiere etape est destinee & decorreler 
I'image en transformant ses pixels en coefficients d'ondelettes repartis dans 
des sous-bandes se situant a plusieurs niveaux de resolution. Une deuxieme 
10 6tape est destinee a quantifier les coefficients d'ondelettes obtenus lors de la 
premiere etape. Une troisieme 6tape est destinee a effectuer un codage 
entropique des coefficients d'ondelettes quantifies lors de la deuxieme 6tape. 

Comme le sait I'homme de Tart, ce type de codage ne peut etre 
optimal qu'a condition que trois hypotheses soient verifiees : la transformee en 
15 ondelettes doit etre orthogonale, les coefficients d'ondelettes doivent etre 
codes par un codage entropique d'ordre 0 dans lequel on n'exploite aucune 
dependance statistique entre coefficients, et on doit se trouver dans les limites 
de la « quantification a haute resolution », dans laquelle les densites de 
probability des coefficients d'ondelettes varient peu a I'interieur d'une boTte de 
20 quantification. 

Lorsque ces hypotheses ne sont pas verifiees, il faut en theorie faire 
une allocation de debit complete, en choisissant un pas de quantification a 
priori different pour chaque sous-bande. Dans la pratique, le codage demeure 
qualitativement valable lorsque les hypotheses ne sont plus formellement 
25 verifiees. C'est notamment le cas lorsque la transformee en ondelettes n'est 
pas orthogonale, mais bi-orthogonale, ou lorsque les coefficients d'ondelettes 
d'une meme sous-bande ne sont pas codes independamment du fait qu'ils ne 
sont pas independants. De meme, lorsque I'hypothese de quantification a 
haute resolution n'est pas verifiee, on peut modifier legerement le codeur en 
30 remplagant son quantificateur uniforme par un quantificateur a intervalle central 
(intervalle de quantification centre en zero ou Zone Morte) de taille double. 
Cela resulte du fait que la distribution des coefficients d'ondelettes dans 
chaque sous-bande est laplacienne (ou gaussienne gen6ralisee) avec un pic 
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en 0 ou elle ne peut plus etre raisonnablement consideree comme quasi 
constante surtout aux debits faibles (c'est-a-dire aux forts taux de 
compression). II est ici rappele qu'une fonction laplacienne est un cas 
particulier d'une gaussienne generalisee. 
5 Pour tenter d'ameliorer le codage dans les conditions precitees 

(transformation bi-orthogonale ou transformation ne conduisant pas a 
I'independance de ses coefficients) il a ete propose de tirer partie de la 
dependance entre coefficients d'ondelettes voisins dans une meme sous- 
bande en realisant un codage entropique contextuel prenant en compte une 
10 densite de probability conditionnelle en lieu et place d'une simple densite de 
probability. Cette densite conditionnelle modelise la loi d'apparition d'un 
coefficient d'ondelettes conditionnellement a son voisinage. 

Dans ce qui suit, on entendra par « dependance intra-bande » une 
dependance entre un coefficient q ( k,k' (ou W J( k[n,m]) et des coefficients voisins 
is de meme echelle ou sous-bande Wjk, comme par exemple q.k-i.k* et q^k'-i ou 
qmv et q,k-i,k-i- De telles dependances intra-bande sont notamment exploitees 
par les codeurs qui effectuent un codage des positions des coefficients 
significatifs en mode « run-length » ainsi que par les codeurs de type 
EBCOT/JPEG2000. 

20 Par ailleurs, dans ce qui suit, on entendra par «dependance inter- 

bandes » une dependance entre un coefficient q^x e t un coefficient place a la 
meme position a Techelle superieure, comme par exemple q+i^m- De telles 
dependances inter-bandes sont notamment exploitees par les codeurs dits 
« en arbre de zeros » (developpes par Shapiro) et le codeur SPIHT (developpe 

25 par Said et Pearlman). 

Ces dependances d6coulent de la representation de chaque 
coefficient d'ondelettes quantifie Q(q,k,k ,n ) sous la forme d'une chaTne 
d'informations comportant son bit de coefficient significatif (N] l k,k ,n = 0 ou 1), et 
lorsque ce bit vaut un (1 ), son signe (Sj,k ( k ,n qui vaut ±1 ) et sa valeur absolue 

30 (A jXk - n ). 

Par exemple, le codage par arbre de zeros exploite la dependance 
statistique entre les bits de coefficient significatif (Nj ik| k ,n )» tandis que le codage 
de type EBCOT exploite la dependance statistique entre les bits de coefficient 
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signiflcatif (N jik ,k ,n ) et la dependance entre signes S jikik . n de coefficients 
d'ondelettes voisins quantifies avec differents pas. 

Ces codages contextuels mettent en ceuvre des methodes de 
traitement particulierement complexes qui imposent aux codeurs d'apprendre 
5 autant d'histogrammes qu'il y a de contextes, ce qui limite notablement leur 
Vitesse d'adaptation et interdit leur utilisation dans le cadre de la compression a 
haut debit d'entree (images de grande taille) et en temps reel, notamment lors 
d'observations spatiales. 

[.'invention a done pour but d'ameliorer la situation. 
!° Elle propose a cet effet un precede dedie a la compression de 

donnees d'images numeriques par transformee en ondelettes et comprenant 
une etape consistant a transformer des donnees d'image numeriques en 
coefficients d'ondelettes repartis dans des sous-bandes, une etape consistant 
a quantifier les coefficients d'ondelettes, et une etape consistant a coder de 
is facon entropique les coefficients d'ondelettes precedemment quantifies. 

Ce precede se caracterise par le fait, d'une part, qu'il comprend, entre 
les etapes de transformation et de quantification, une etape consistant a 
estimer, pour chaque coefficient de chaque sous-bande d'une image, de 
premier et second jeux de parametres de prediction, associes respectivement 
20 a des directions dites « nord-sud » et « ouest-est », en fonction des valeurs des 
coefficients d'ondelettes de ses voisins nord ou ouest, et d'autre part, qu'a 
I'etape de codage entropique, d'une part, on determine pour chaque coefficient 
de chaque sous-bande de I'image des valeurs de prediction de I'esperance et 
de la largeur d'une fonction laplacienne (ou gaussienne generalisee), 
25 representative de sa densite de probability, en fonction des premier et/ou 
second jeux de parametres de prediction et du coefficient d'ondelettes quantifie 
de son voisin nord ou de son voisin ouest, et d'autre part, on code de fagon 
entropique les coefficients d'ondelettes quantifies a I'aide des esperances et 
largeurs associees determinees. 
so Grace aux estimations des parametres de prediction, la phase de 

prediction de I'etape de codage entropique est notablement simplifies, ce qui 
penmet de compresser (ou comprimer) des donnees d'image numeriques a 
haut debit et en temps reel, notamment lors d'observations spatiales. 
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Le procede selon I'invention peut comporter d'autres caracteristiques 
qui pourront etre prises separement ou en combinaison, et notamment : 

- lorsque cela s'avere possible, on peut determiner lors de I'etape de codage 
entropique, pour la fonction laplacienne de chaque coefficient d'ondelettes 
quantifie, des valeurs de prediction d'esperance nord et ouest et des valeurs 
de prediction de largeur nord et ouest. La valeur de prediction d'esperance 
nord (ou ouest) est alors preferentiellement definie par le produit entre un 
premier parametre de prediction du premier (ou second) jeu et la valeur du 
coefficient d'ondelettes quantifie voisin nord (ou ouest), et la valeur de 
prediction de largeur nord (ou ouest) est definie par la somme d'un 
deuxieme parametre de prediction du premier (ou second) jeu et du produit 
d'un troisieme parametre de prediction du premier (ou second) jeu et de la 
valeur absolue du coefficient d'ondelettes quantifie voisin nord (ou ouest). 
Dans ce cas, quatre cas doivent §tre envisages selon I'environnement du 
coefficient d'ondelettes concerne. En I'absence de coefficient voisin nord, on 
identifie les valeurs de prediction d'esperance et de largeur de la fonction 
laplacienne aux valeurs de prediction d'esperance ouest et de largeur ouest. 
En I'absence de coefficient voisin ouest, on identifie les valeurs de prediction 
d'esperance et de largeur de la fonction laplacienne aux valeurs de 
prediction d'esperance nord et de largeur nord. En I'absence de coefficients 
voisins ouest et nord, on effectue un codage laplacien fixe. Enfin, en 
presence de coefficients voisins ouest et nord, on identifie les valeurs de 
prediction d'esperance et de largeur de la fonction laplacienne 
respectivement a la demi-somme des valeurs de prediction d'esperance 
nord et ouest et a la demi-somme des valeurs de prediction de largeur nord 
et ouest, 

- on peut estimer les parametres de prediction des premier et second jeux par 
regression. 

• Par exemple, on peut estimer les premiers parametres de prediction des 
premier et second jeux par une premiere regression obtenue en 
resolvant une equation de type W[n+1,m] = 0.W[n,m], ou W[n,m] est la 
valeur d'un coefficient d'ondelettes et W[n+1,m] est la valeur du 
coefficient d'ondelettes voisin de W[n,m], et 0 est un premier parametre 
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de prediction. On peut par exemple resoudre la premiere equation a 
I'aide de la technique des moindres carres, le premier parametre de 
prediction d'un premier ou second jeu etant alors egal au rapport entre 
les covariance et variance non normalisees des valeurs des coefficients 
d'ondelettes associes au premier ou second jeu. 
• De meme, on peut estimer, d'une part, les deuxieme et troisieme 
parametres de prediction du premier jeu par une seconde regression 
obtenue en resolvant un premier systeme d'equations de type |W[n,m] 
- 0 N .W[n,m-1] | = a N + fin |W[n,m-1] |, ou W[n,m] est la valeur d'un 
coefficient d'ondelettes, W n ,m.i est la valeur du coefficient d'ondelettes 
voisin nord de \N n , m , et 0 N , a N et fin sont respectivement des premier, 
deuxieme et troisieme parametres de prediction du premier jeu, et 
d'autre part, les deuxieme et troisieme parametres de prediction du 
second jeu par une seconde regression obtenue en resolvant un second 
systeme d'equations de type |W[n,m] - 0 o .W[n-1,m] | = a 0 + fi 0 |W[n- 
1,m]|,o0 W n -i. m est la valeur du coefficient d'ondelettes voisin ouest de 
W n ,m, et 0 O , ao et fi 0 sont respectivement des premier, deuxieme et 
troisieme parametres de prediction du second jeu. On peut par exemple 
resoudre les premier et second systemes d'equation par la technique 
des moindres carres, 

- lorsqu'un coefficient d'ondelettes quantifie voisin presente une valeur nulle, 
on peut le substituer par le produit du pas de quantification, utilise lors de 
I'etape de quantification, et d'une constante choisie, strictement superieure a 
zero et inferieure ou egale a un (1), et plus preferentiellement egale a 1/3, 

- lors de i'etape de codage entropique, on peut proceder a la discretisation de 
la fonction laplacienne apres avoir determine les valeurs de prediction de 
I'esperance et de la largeur. En variante, on peut proceder au codage 
entropique par iteration apres avoir determine les valeurs de prediction de 
I'esperance et de la largeur. 

L'invention porte egalement sur un dispositif dedie a la compression de 
donnees d'images numeriques par transformee en ondelettes et comprenant 
des moyens de traitement munis, d'une premiere part, d'un module charge de 
transformer les donnees d'image numeriques qu'il recoit en coefficients 
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d'ondelettes repartis dans des sous-bandes, d'une deuxieme part, d'un module 
charge de quantifier les coefficients d'ondelettes, et d'une troisieme part, d'un 
module charge de coder de fagon entropique les coefficients d'ondelettes 
quantifies. 

5 Ce dispositif se caracterise par le fait, d'une part, que ses moyens de 

traitement comprennent egalement un module charge d'estimer, pour chaque 
coefficient de chaque sous-bande de I'image, de premier et second jeux de 
parametres de prediction, associes respectivement a des directions dites 
« nord-sud » et « ouest-est », en fonction des valeurs des coefficients 
10 d'ondelettes de ses voisins nord ou ouest, et d'autre part, que son module de 
codage entropique est charge, d'une part, de determiner, pour chaque 
coefficient de chaque sous-bande de I'image, des valeurs de prediction de 
I'esperance et de la largeur d'une fonction laplacienne (ou gaussienne 
generalisee), representative de sa densite de probability en fonction des 
15 premier et/ou second jeux de parametres de prediction et du coefficient 
d'ondelettes quantifie de son voisin nord ou de son voisin ouest, et d'autre part, 
de coder de fagon entropique les coefficients d'ondelettes quantifies a I'aide 
des esperances et largeurs associees determinees. 

Le dispositif selon I'invention peut comporter d'autres caracteristiques 
20 qui pourront etre prises separement ou en combinaison, et notamment : 

- son module de codage entropique peut etre agence de maniere a 
determiner pour la fonction laplacienne de chaque coefficient d'ondelettes 
quantifie, lorsque cela s'avere possible, des valeurs de prediction 
d'esperance nord et ouest et des valeurs de prediction de largeur nord et 
25 ouest. La valeur de prediction d'esperance nord (ou ouest) est alors 
preferentiellement definie par le produit entre un premier parametre de 
prediction du premier (ou second) jeu et la valeur du coefficient d'ondelettes 
quantifie voisin nord (ou ouest), et la valeur de prediction de largeur nord (ou 
ouest) est definie par la somme d'un deuxieme parametre de prediction du 
30 premier (ou second) jeu et du produit d'un troisieme parametre de prediction 
du premier (ou second) jeu et de la valeur absolue du coefficient 
d'ondelettes quantifie voisin nord (ou ouest). Dans ce cas, le module de 
codage entropique est preferentiellement agence pour traiter quatre 



2005/008886 PCT/FR2004/001777 

7 

situations distinctes. En I'absence de coefficient voisin nord, il identifie les 
valeurs de prediction d'esperance et de largeur de !a fonction lapiacienne 
aux valeurs de prediction d'esperance ouest et de largeur ouest. En 
I'absence de coefficient voisin ouest, il identifie les valeurs de prediction 
d'esperance et de largeur de la fonction lapiacienne aux valeurs de 
prediction d'esperance nord et de largeur nord. En I'absence de coefficients 
voisins ouest et nord, il effectue un codage laplacien fixe. Enfin, en presence 
de coefficients voisins ouest et nord, il identifie les valeurs de prediction 
d'esperance et de largeur de la fonction lapiacienne respectivement a la 
demi-somme des valeurs de prediction d'esperance nord et ouest et a la 
demi-somme des valeurs de prediction de largeur nord et ouest, 
son module d'estimation peut etre agence de maniere a estimer les 
parametres de prediction des premier et second jeux par regression. 

• Par exemple, le module d'estimation peut etre charge d'estimer les 
premiers parametres de prediction des premier et second jeux a I'aide 
d'une premiere regression obtenue en resolvant une premiere equation 
de type W[n+1,m] = 0.W[n,m], ou W[n,m] est la valeur d'un coefficient 
d'ondelettes et W[n+1,m] est la valeur du coefficient d'ondelettes voisin 
de W[n,m], et 6 est un premier parametre de prediction. Dans ce cas, le 
module d'estimation peut par exemple resoudre la premiere equation par 
la technique des moindres carres. Le premier parametre de prediction 
d'un premier ou second jeu etant alors egal au rapport entre les 
covariance et variance non normalises des valeurs des coefficients 
d'ondelettes associes au premier ou second jeu. 

• Par ailleurs, le module d'estimation peut etre charge, d'une part, d'estimer 
les deuxieme et troisieme parametres de prediction du premier jeu a 
I'aide d'une seconde regression obtenue en resolvant un premier 
systeme d'equations de type | W[n,m] - 0 N .W[n,m-1] | = a N + fa | W[n,m- 
1] | , ou W[n,m] est la valeur d'un coefficient d'ondelettes, W n ,m-i est la 
valeur du coefficient d'ondelettes voisin nord de W n , m , et 0 N) a N et fa sont 
respectivement des premier, deuxieme et troisieme parametres de 
prediction du premier jeu, et d'autre part, d'estimer les deuxieme et 
troisieme parametres de prediction du second jeu a I'aide d'une seconde 
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regression obtenue en resolvant un second systeme d'equations de type 
I W[n,m] - 0 o .W[n-1 ,m] | = a 0 + fa | W[n-1 ,m] | , ou W n . 1>m est la valeur du 
coefficient d'ondelettes voisin ouest de W n>m , et e 0 , a 0 et p Q sont 
respectivement des premier, deuxieme et troisieme parametres de 
5 prediction du second jeu. Dans ce cas, le module d'estimation peut par 

exemple resoudre les premier et second systemes d'equation a I'aide de 
la technique des moindres carres, 

- son module de codage entropique peut etre agence, lorsqu'un coefficient 
d'ondelettes quantifie voisin presente une valeur nulle, de maniere a le 

10 substituer par le resultat d'un produit entre le pas de quantification, utilise par 
le module de quantification, et une constante choisie, strictement superieure 
a zero et inferieure ou egale a un (1 ), et preferentiellement egale a 1/3, 

- son module de codage entropique peut etre agence de maniere a proceder 
a une discretisation de la fonction laplacienne apres avoir effectue les 

15 predictions d'esperance et de largeur. En variante, le module de codage 
entropique peut etre agence de maniere a proceder au codage par iteration 
apres avoir effectue les predictions d'esperance et de largeur. 

L'invention porte egalement sur un compresseur/decompresseur de 
donnees d'image numeriques, ou codec, comprenant un dispositif de 
2 0 compression du type de celui presente ci-avant. 

D'autres caracteristiques et avantages de l'invention apparaTtront a 
I'examen de la description detaillee ci-apres, et des dessins annexes, sur 
lesquels : 

- la figure 1 illustre schematiquement, sous la forme d'un diagramme bloc 
25 fonctionnel, un dispositif de compression de donnees d'image numeriques 

selon l'invention, et 

- la figure 2 illustre une decomposition d'image en sous-bandes de niveaux 
de resolution differents. 

Les dessins annexes pourront non seulement servir a completer 
30 l'invention, mais aussi contribuer a sa definition, le cas echeant. 

L'invention a pour objet de permettre la compression et la 
decompression de donnees d'images numeriques par transformee en 
ondelettes. 
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La compression (ou le codage) par transformee en ondeiettes consiste 
classiquement a decomposer une image ! en sous-bandes Wj, k dites LH, HL, 
HH et LL, puis a transformer les donnees d'image numeriques l[n,m], qui 
definissent I'intensit6 et la couleur des pixels de I'image I, en ondeiettes, et plus 

5 precisement en un tableau [n,m] de coefficients d'ondelettes, puis a quantifier 
les coefficients d'ondelettes avant de les coder de facon entropique. 

n et m sont ici des indices entiers qui varient dans les limites de 
I'image : n designe la colonne et appartient a I'intervalle [0, n max - 1] et m 
designe la ligne et appartient a I'intervalle [0, ow - 1]. Une sous-bande W jik est 

10 done definie par un tableau bidimensionnel de coefficients d'ondelettes 
Wj,k[n,m] ou les indices n et m varient dans des sous-intervalles [0, nmax / 2 j - 1] 
et [0, m max / 2 J - 1]. On note alors N = n^x . m max le nombre de pixels de I'image 
I, et N jik = N / 2 2j le nombre de coefficients d'ondelettes dans une sous-bande 
W jlk . 

is Les sous-bandes sont codees separement les unes des autres, la 

sous-bande passe bas etant codee avec un codage differentiel de type DPCM 
et les differences resultantes pouvant etre codees avec une methode similaire 
a celle utilisee pour une sous-bande passe-haut. 

Par ailleurs, on designe par Wj, k [n,m-1] et W j>k [n-1,m] les voisins dits 

20 «nord» et «ouest» d'un coefficient d'ondelettes W j>k [n,m]. Les coefficients 
d'ondelettes a coder sont parcourus dans un ordre particulier (ordre 
lexicographique croissant sur le couple (n,m)) afin de garantir qu'un coefficient 
Wj, k [n,m] sera toujours code apres ses voisins nord et ouest qui servent a la 
determination de sa valeur. 

25 La compression selon I'invention differe d'une compression par 

transformee en ondeiettes classique notamment par le fait qu'elle comprend 
entre son etape de transformation (classique) et son etape de quantification 
(classique), une etape dans laquelle on estime, pour chaque coefficient 
W jik [n,m] de chaque sous-bande W jik d'une image I, de premier et second jeux 

30 de parametres de prediction, associes respectivement a des directions nord- 
sud et ouest-est, en fonction des valeurs des coefficients d'ondelettes de ses 
voisins nord W jik [n,m-1] ou ouest W jik [n-1 ,m]. 

Ces premier et second jeux de parametres de prediction sont ensuite 
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utilises lors de I'etape de codage entropique pour determiner pour chaque 
coefficient d'ondelettes Wj^n.m] de chaque sous-bande W jik de 1'image ! des 
valeurs de prediction de I'esperance (ou moyenne) jj et de la largeur a d'une 
fonction (ou distribution) de type laplacienne representative de sa densite de 
5 probability en fonction d'un coefficient d'ondelettes quantifie voisin (nord ou 
ouest) : 

2o-(c') 

ou c represente un coefficient d'ondelettes Wj, k [n,m], c' represente soit le 
coefficient d'ondelettes voisin nord Wj ik [n,m-1] de c, soit le coefficient 
10 d'ondelettes voisin ouest vV jik [n-1 ,m] de c, et /y(c') et ct(c') represented 
respectivement I'esperance et la largeur en fonction du coefficient d'ondelettes 
voisin c' considere. 

II est ici rappele que I'esperance u et la largeur a d'une distribution 
laplacienne de densite P sont definies par les relations suivantes : 

15 n = jxP(x)dx 

cr= jjjc - [i\P(x)dx 

On va maintenant decrire plus en detail les differentes etapes de la 
compression selon invention. 

Comme indique ci-avant, la compression commence par une etape de 
20 transformation des pixels de I'image I en coefficients d'ondelettes W j)k [n,m]. 
Cette etape de transformation etant tout a fait classique, elie ne sera pas 
decrite ici. Pour la mettre en ceuvre on peut par exemple utiliser la transformee 
en ondelettes bi-orthogonales 7-9 de Daubechies/Antonini. Cette transformee 
n'est pas rigoureusement orthogonale, mais elle I'est quasiment, ce qui permet 
25 de conserver I'energie des images a quelques pour cent pres. 

A la fin de cette etape de transformation, on dispose done pour une 
sous-bande donnee W j>k de I'ensemble de ses coefficients d'ondelettes 
W jik [n,m]. 

L'etape de transformation est suivie d'une etape d'estimation de 
30 parametres de prediction pour chaque coefficient d'ondelettes W j)k [n,m] de 
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chaque sous-bande \%. Dans ce qui suit, et afin de simplifier les notations on 
ne precisera plus les indices j et k des sous-bandes. Par consequent, on 
designera par W[n,m] un coefficient d'ondelettes appartenant a une sous- 
bande W jtk [n,m], qu'ii s'agisse d'une sous-bande LH, HH, ou HL, ou encore LL 
5 (par exemple codee en DPCM par un parcours en boustrophedon). 

Comme indique precedemment, I'etape d'estimation est destinee a 
determiner des premier et second jeux de parametres de prediction qui 
permettront ulterieurement (lors de I'etape de codage entropique) de predire 
les valeurs de I'esperance (ou moyenne) u et de la largeur a d'une fonction (ou 
10 distribution) laplacienne representative de la densite de probability d'un 
coefficient d'ondelettes quantifie compte tenu de la valeur prise par le 
coefficient d'ondelettes quantifie de I'un au moins de ses voisins nord et ouest. 

L'esperance fj et la largeur a sont definies par rapport aux parametres 
de prediction des premier et second jeux. Par exemple, chaque jeu comporte 
is trois parametres de prediction 9 (premier), a (deuxieme) et J3 (troisieme), et les 
relations liant ces parametres aux esperance jj et largeur a prennent les 
formes suivantes : 

u = 0.C (1) 
a = a + J3.c' (2) 

20 ou, comme indique precedemment c' designe le contexte representant soit le 
coefficient d'ondelettes voisin nord W[n,m-1] de c (ou W[n,m]), soit le coefficient 
d'ondelettes voisin ouest W[n-1 ,m] de c. 

En fait, comme on le verra plus loin, on utilise deux predictions 
d'esperance et de largeur selon le voisin considere (nord ou est) : une 
25 prediction nord porte sur I'esperance nord un et la largeur nord ct n , et une 
prediction ouest porte sur I'esperance ouest u N et la largeur ouest a N . Le 
premier jeu est alors associe a la prediction nord et comporte les parametres 
0 N , o N et /?n, tandis que le second jeu est associe a la prediction ouest et 
comporte les parametres 6 0 , a Q et p 0 . 

30 Les relations (1) et (2) presentees ci-dessus, a titre d'exemple, 

definissent des regies de calcul qui permettent d'estimer simplement les 
parametres 9,aet0 (nord et ouest) avec des regressions. 

Preferentiellement, on commence par estimer 0 N et 6 0 par une 
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premiere regression. Celle-ci s'obtient en resolvant une equation du type : 

Cn +1 = e.C 

En utilisant, par exemple, la technique des moindres carres, on obtient 
alors le premier parametre 6 pour chaque jeu (nord ou ouest) : 

5 0 = -s- 



n 

La relation ci-dessus peut etre reecrite comme suit, pour les jeux nord 
et ouest : 

N Var 

00=^- 
0 Var 

10 ou Covn, Cov 0 et Var sont respectivement les covariances Nord et Ouest et la 
variance, non normalisees, donnees par les relations suivantes : 

Var = Y J W[n,m] 2 

Cov N = £ W[n, m}v[n, m - 1] 

Cov o = ^W[n,m}v[n - l,m] 

is Puis, on procede a Pestimation des parametres de prediction nord an 

et fin par une seconde regression. Celle-ci s'obtient en r§so!vant un premier 

systeme d'equations du type : 

| W[n,m] - 0 N .W[n,m-1] | = a N + 0n | W[n,m-1] | (3) 

En imposant alors les definitions suivantes : 

20 M n =£l (4) 

«,/» 

Mn=M 2X ^\w{n,m-\} (5) 

M 22 =Z\ w i n > m - 1 Y (6) 

Y\ = jyM- e N W[n,m - 1] (7) 

^Zrtv»]-^h-iM".«i-i] (8) 



WO 2005/008886 PCT/FR2004/001777 

13 

on obtient les parametres de prediction an et 0 H en resolvant le systeme : 



M 2\ M 22 j 



Les parametres de prediction ouest a 0 etfio sont calcules de la meme 
maniere, en resolvant un second systeme d'equations du type : 
5 I W[n,m] - 0 o .W[n-1 ,m] | = a Q + p 0 | W[n-1 ,m] | (1 0) 

Cela revient a resoudre le systeme defini par la relation (9), en utilisant 
des coefficients definis par les relations (4) a (8) dans lesquelles 0 N est 
remplace par do et W[n,m-1] est remplace par W[n-1,m], et en remplacant a N 
par a 0 et^ N par/? 0 . 

10 Une fois que I'on a determine les parametres de prediction nord (0 N , a u 

et /?n) et ouest (6 0 , a 0 et fio) de chaque coefficient d'ondelettes W j|k [n,m] d'une 
sous-bande Wj, k) on passe a I'etape de quantification. 

Cette etape de quantification etant tout a fait classique, elle ne sera 
pas decrite ici. Pour la mettre en oeuvre on peut par exemple utiliser un 



W[n,m] 



oD T est la 



15 quantificateur scalaire standard du type Q(W[n,rn]) = 

largeur du pas de quantification. 

II est simplement rappele que I'etape de quantification consiste a 
quantifier les coefficients d'ondelettes W jik [n,m] de chaque sous-bande W jik a 
partir d'un pas de quantification T soit fixe, soit determine par un module de 
20 planification. 

En raison de la nature meme des precedes de compression, le debit 
de sortie (du dispositif de compression) est directement proportionnel a 
I'entropie du signal d'entree (donnees d'images numeriques a compresser). Or, 
la quantite d'information a transmettre varie generalement au cours du temps 

25 puisque les caracteristiques du signal d'entree varient elles-memes avec le 
temps. Par consequent, dans le but d'optimiser ('utilisation du canal de 
transmission, on minimise les fluctuations de debit en sortie du dispositif de 
compression en effectuant une regulation de debit. Cela est meme 
indispensable lorsque le canal de transmission necessite un debit constant. 

30 Pour permettre une telle regulation de debit en temps reel, il est 

imperatif de mettre en oeuvre des algorithmes de planification de debit rapides 
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(exempts d'optimisation sous contraintes et exempts alterations) de maniere a 
ne pas ralentir I'acquisition des images lors du defilement du satellite ou ; dans 
le cadre du codage video, pour ne pas generer d'effet de saccade. 

Deux types de planification de debit ont ete proposes a cet effet : la 
s planification avec allocation de debit et la planification de debit pure. 

La planification avec allocation de debit consiste a laisser le module de 
planification, d'une part, choisir independamment les pas de quantification des 
differentes sous-bandes de sorte que le debit total soit egal a une consigne de 
debit, et d'autre part, choisir la facon dont ce debit doit etre reparti dans les 
10 differentes sous-bandes. 

La planification de debit pure consiste a Her les pas de quantification 
des differentes sous-bandes (en general ils sont identiques). Dans ce cas, le 
module de planification se contente de choisir le pas de quantification unique 
pour toutes les sous-bandes qui permettra d'atteindre une consigne de debit 
is fixee. 

Par ailleurs, afin de relacher la contrainte de debit fixe on peut 
controler la qualite de la compression (et done de la decompression) pour 
imposer une Constance de la qualite des images en sortie du dispositif de 
compression au gre du defilement du satellite ou de I'arrivee des images en 
20 codage video. 

A cet effet, il a ete propose un procede de planification de qualite 
consistant a determiner un pas de quantification uniforme T permettant 
d'obtenir une qualite de codage donnee (generalement fixee en dB). Ce 
procede consiste tout d'abord a calculer la distorsion totale D correspondant a 
un pas de quantification donne T, puis a determiner par dichotomie le pas de 
quantification T permettant d'atteindre une consigne de distorsion donnee D. 

D'une facon plus precise, dans le domaine de I'allocation de debit ou 
du contrSle de qualite il existe principalement trois types de procede. 

Un premier type concerne les procedes qui realisent une optimisation 
sous contraintes par une minimisation d'un lagrangien sous la forme J = D+ AR 
(ou D est la distorsion et R le d6bit), avec pour contrainte R = Ro (ou Ro est un 
debit de consigne). 

Un deuxteme type concerne les procedes fondes sur une moderation 



25 
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de la chame de compression pour etablir une regie qui lie le debit R au pas de 
quantification global T ou bien la distorsion D au debit R. 

Ces deux precedes ne peuvent cependant pas etre utilises lorsque Ton 
a besoin de travailler avec des decompositions non uniforme, comme par 
5 exemple lors d'une re-decomposition de sous bandes pour la prise en compte 
de la FTM (Fonction de Transfert de Modulation) en vue d'une deconvolution 
ou d'une allocation par sous-bandes dans le cas d'images multicomposantes. 

Un troisieme type concerne les precedes fondes sur une modelisation 
de la chaTne de compression pour etablir une regie qui lie le debit R au pas de 
10 quantification local T a chacune des sous-bandes ou bien la distorsion D au 
debit R. Un precede de ce type, adapte a I'allocation de debit par sous-bandes 
a I'aide d'un codeur arithmetique d'ordre 1, est notamment decrit dans le 
document FR 02/05724. 

Un autre but de I'invention est d'ameliorer le resultat procure par le 
is compresseur (ou codeur) decrit dans le document precite en etablissant une 
regie optimale de planification du debit ou de la qualite d'image d'un 
compresseur par transformee en ondelettes grace au gain de codage 
engendre par la prise en compte des dependances entre coefficients 
d'ondelettes intra-bande. 
20 Pour atteindre ce but, il est tout d'abord necessaire de modeliser le 

comportement d'une sous-bande en termes de debit R et de distorsion D, en 
fonction des parametres de prediction estimes 0, a et de la fonction 
laplacienne conditionnelle, determines lors de I'etape d'estimation presentee ci- 
avant. 

25 Pour effectuer cette modelisation, il est possible d'utiliser 

simultanement, pour chaque sous-bande, les deux directions Nord-Sud et 
Ouest-Est, ou en d'autres termes les premier (0 Nl <r N , /? N ) e t second (0 O> a Q , fo) 
jeux de parametres de prediction. Mais, cela imposerait d'effectuer la carte 
d'une fonction de grande dimension, et de plus, on peut montrer que le codage 

30 d'un coefficient d'ondelettes connaissant un seul coefficient d'ondelettes non 
quantifie voisin (nord ou ouest) donne un resultat de codage semblable a un 
codage utilisant les coefficients d'ondelettes non quantifies des voisins nord et 
ouest. Par consequent, dans ce qui suit on se place dans une situation dans 
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laquelle on n'utilise qu'un seul coefficient d'ondelettes non quantifie voisin (nord 
ou ouest). 

Dans cette situation « monodirectionnelle », on peut definir une chaTne 
de dependance, dite « de Markov », par la relation suivante : 

S P(X n |X„- 1 ) =- 1 re H^-^.-.|/(^.-,|) 

On peut alors montrer que si pour un jeu de parametres de prediction 
donnes 9, a et p, I'esperance et la variance de la variable aleatoire X n restent 
bornees quand n croTt, alors la suite de variables aleatoires est ergodique. 

Pour calculer le debit R et la distorsion D d'une telle suite, pour un pas 
10 de quantification T donne, on peut utiliser une methode de Monte-Carlo 
consistant a faire des tirages d'une sequence de la chaTne de Markov, puis a 
estimer les esperances de debit R et de distorsion D par symbole en effectuant 
une moyenne sur des termes consecutifs de la sequence. 

Ces calculs de debit R et de distorsion D presentent un sens lorsque la 
is chaine de Markov est ergodique, et notamment des que sa largeur (ou 
variance) est bornee. 

Les tirages de la chaTne de Markov peuvent etre faits de la maniere 
suivante. On calcule tout d'abord par recurrence Xn + i en fonction de X n par une 
formule du type : 
2 o X n+ i = S n (or + ^ | X n | ) logUn + 0X n 

ou S n et U n sont des variables aleatoires independantes, chaque S n 
valant +1 et -1 avec des probabilites respectives 1/2 et 1/2, et U n ayant une 
densite uniforme sur I'intervalle [0,1]. 

On doit ensuite estimer des courbes de debit et de distorsion. Pour ce 
25 faire, on doit tout d'abord estimer les relations qui lient le pas de quantification 
T au debit R et a la distorsion D, pour une gamme de parametres de prediction 
9,ae\p pertinente. 

Par definition, la gamme des valeurs pertinentes est celle pour laquelle 
la chaTne de Markov est ergodique. On peut remarquer que I'ergodicite ne 
30 depend pas du parametre de prediction a qui est un parametre d'echelle. En 
effet, si une chaTne de variables aleatoires X n est une chaTne de Markov pour 
les parametres de prediction 6, a et p, alors la chaTne AX n est une chaTne de 
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Markov pour les parametres 9, Act et 0. On peut done se limiter au cas a = 1 . 

Par ailleurs, i! est naturel de limiter les vaieurs de 0 a I'intervalle [0,1], 
puisque 0 ne peut pas etre negatif. Si 0 est plus grand que 1, il est 
vraisemblable que la variance de la distribution associee tende vers I'infini (par 

s un etalement progressif). De meme, 9 doit presenter une valeur absolue 
inferieure a 1, sans quoi le processus risquerait d'avoir un comportement de 
marche aleatoire, lequel n'est pas ergodique. En revanche, le domaine des 
parametres de prediction 9 licites (les autres parametres de prediction etant 
fixes) est necessairement symetrique, puisque si une suite X„ est de 

10 parametres 9, a et 0, la suite (-1) n X n est de parametres -9, a et 0. Par un 
raisonnement analogue, on peut meme montrer que le domaine en question 
est un intervalle de la forme ]-9 0 , 9 0 [. 

Faute de pouvoir definir clairement par des methodes analytiques un 
domaine des couples (/?, 9) pour lequel la chame de Markov est ergodique, on 

is en definit un sous-ensemble pour lequel la largeur estimee de la variable 
aleatoire E( | Xn | ) est bornee par une limite fixee. 

Lorsque I'on trace la carte de la variable aleatoire E(X n ) en fonction des 
parametres de prediction 0 et 9, on constate que I'estimation d'une esperance 
devient de plus en plus couteuse (en nombre d'iterations) au fur et a mesure 

20 que I'on s'approche de la non-ergodicite. Par ailleurs, en tracant des sections 
de cette fonction de deux variables et de son inverse, on peut constater que 
son inverse est sensiblement lineaire en 0 : 



ou A et B sont des fonctions a vaieurs positives. 
25 Par consequent, pour chaque section de cette fonction a deux 

variables, a 9 constant, la pente et I'ordonnee a I'origine peuvent etre estimees 
par une regression tres stable. Par exemple, les fonctions A(0) et B(0) peuvent 
etre parametrees par decomposition dans une base de Deslauriers-Dubuc. II 
est done globalement possible d'estimer les limites d'ergodicite avec une 
30 grande precision. 

Le parametrage des fonctions A(9) et B(0) au voisinage de un (1) sont 
plus precis si on decompose dans une base de Deslauriers-Dubuc non pas les 
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fonctions A(0) et B(e), mais plut6t A(*) 2 et B(ef. En effet les graphes de ces 
fonctions semblent montrer qu'elles ont un zero de la forme (1 - 9) V2 en un (1) 
pour A, et pour une valeur 9 > 1 pour B. La decomposition dans la base de 
Deslauriers-Dubuc de chaque fonction coQte neuf coefficients, 
s Si I'on appelle R(0, a, 0, T) la fonction de debit d'une sous-bande en 

fonction des parametres de prediction 9, a et 0 et du pas de quantification T, on 
peut la definir par la relation suivante : 

R(9, a, 0, T) = E(-P(Q(X n ) | Xm) log 2 P(Q(X n ) | X^)) 
ou X n est une chaTne de Markov laplacienne de parametres 9, a et 0, 
io et Q un quantificateur de pas T. 

Le debit R peut etre estime par des moyennes de Monte-Carlo 
calculees sur des tirages de la chaTne de Markov, done par moyennage de la 
suite : 

Rn = -PnlOg 2 P n 



15 OU 



Ces moyennes peuvent §tre estimees pour des parametres de 
prediction 6 et 0 qui sont dans le domaine d'ergodicite identifies plus haut. On 
peut utiliser I'homogeneite de la fonction de debit R pour reduire la 
2 o dimensionnalite du probleme, du fait que R(0, a, 0, T) = R(0, 1 , 0, ji a ). 

Si I'on appelle maintenant D(9, a, 0, T) la fonction de distorsion de la 
chaTne de Markov, qui servira a modeliser la distorsion D d'une sous-bande 
quantifiee en fonction des parametres de prediction 9, a et 0 et du pas de 
quantification T, on peut la definir par la relation suivante : 
25 D(0,a,/?,T)=E(|Q(X n )-X n | 2 ) 

ou X n est de nouveau la chaTne de Markov laplacienne de parametres 
9, a et/?, et Q un quantificateur de pas T. 

Ici encore, la distorsion D peut etre estimee par des moyennes de 
Monte-Carlo de la suite D n = |X n - Q(X n )|. II faut en effet connaTtre la 
so distribution stationnaire associee a la chaTne de Markov. Or, il n'en existe pas 
d'expression explicite, de sorte qu'il faut recourir a une methode de Monte- 
Carlo. 
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lei encore, I'homogeneite permet de reduire la dimension du probleme 
D(9, a 1 /3,T) = a z D(9^,/3,T/a) 

Etant donne que ie domaine d'ergodicite contient seulement trois 
5 points singuliers (ou coins), on peut utiliser un maillage triangulaire du domaine 
d'ergodicite concerne et utiliser des fonctions splines triangulaires pour 
representer les fonctions de distorsion D et de debit R sur ce domaine 
d'ergodicite. Cependant, on peut montrer que le coin superieur du domaine est 
un point ou les fonctions D et R peuvent avoir de grandes variations, et qu'il est 
10 preferable d'utiliser un maillage plus fin autour de ce point, comme par exemple 
un maillage carre projete dans le domaine d'ergodicite. 

Pour 9, 0 dans le domaine d'ergodicite maille, a = 1 et la variable T 
(pas de quantification) variant entre 1CJ 3 et 1000, des estimations de la 
distorsion D et du debit R peuvent etre faites par la methode de Monte-Carlo 
is sur des grandes series (typiquement quelques millions de tirages). On obtient 
ainsi une carte de D et R sur un domaine d'ergodicite qui couvre une gamme 
de parametres tres large. 

Ensuite, la fonction inverse (qui associe un pas de quantification T a un 
debit R) peut etre calculee et decomposee dans une base de Deslauriers- 
20 Dubuc, par inversion de I'application 9, fi, T -* R en 9, p, R -» T. 

Dans cette inversion, il faut prendre un parametrage en R soit lineaire, 
soit iogarithmique, afin de representer fidelement la fonction T(R) dans des 
regions ou elle varie tres brutalement. 

On peut ensuite effectuer la prediction des relations entre le debit R, le 
25 pas de quantification T et la distorsion D pour une sous-bande donnee. 

Comme on I'a vu precedemment, pour chaque sous-bande la 
prediction d'un coefficient d'ondelettes se fait preferentiellement a I'aide de 
deux coefficients d'ondelettes voisins (nord et ouest), alors que les modeles de 
debit R et de distorsion D presentes ci-avant n'utilisent qu'un contexte a un seul 
3 o coefficient d'ondelettes voisin. 

Pour associer un modele monodimensionnel a une sous-bande, il faut 
alors reduire la paire de premier et second jeux .de parametres de prediction, 
associes a deux fonctions laplaciennes marginales, a un unique jeu. La 
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reduction utilisee pour la planification de debit est par exemple la suivante : 

a - Ja N a Q 

p 2J2 

2 

5 On peut voir cette reduction comme un moyennage soit arithmetique, 

soit geometrique. Seul & est en plus encore contracts d'un facteur *J2, cense 
modeliser la reduction d'incertitude liee a I'usage d'un contexte double. 

Le choix du facteur *fi peut etre justifie par une approche bayesienne 
simple : deux mesures, bruitees par des bruits independants combines, 
10 peuvent etre combinees (par moyennage) en reduisant I'ecart-type du bruit de 
mesure d'un facteur ^J2. Le choix de cette convention est certes arbitraire, mais 
il permet d'obtenir des resultats raisonnables de prediction pour des images 
classiques. Bien entendu, d'autres choix peuvent etre envisages pour affiner 
cette reduction. 

is Compte tenu de ce qui a ete precise ci-avant, il est desormais possible 

d'elaborer simplement des algorithmes de planification. 

Ainsi, dans le cas d'une planification de debit pure, consistant a 
determiner un pas de quantification uniforme T permettant d'atteindre un debit 
de consigne global pour une image, il suffit de calculer pour chaque sous- 

20 bande une prediction de debit en fonction d'un pas de quantification T, soit 
done un debit total R. Pour ce faire, on determine par dichotomie la valeur du 
pas de quantification T qui rend le debit R, correspondant a T, egal a la 
consigne de debit. 

Dans le cas d'une planification par allocation de d6bit, on cherche a 
25 determiner librement un pas de quantification T pour chaque sous-bande, afin 
de minimiser la distorsion D pour une consigne de debit globale R. Pour ce 
faire, on decoupe le debit total R en trongons de petite taille, et on distribue ces 
troncons de fagon iterative entre les differentes sous-bande, en donnant a 
chaque fois un troncon a la sous-bande qui reduit le plus la distorsion globale 
30 D. Comme indique precedemment, un exemple de mise en ceuvre de cette 
methode est decrit dans le document FR 02/05 724. 
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Dans le cas d'une planification de qualite, on cherche a determiner un 
pas de quantification uniforme T permettant d'obtenir une qualite de codage 
donnee (fixee en dB). Pour ce faire, on calcule la distorsion totale D qui 
correspond a un pas de quantification donne T, puis on determine par 
s dichotomie le pas de quantification T qui permet d'atteindre la consigne de 
distorsion donnee D. 

A la fin de I'etape de quantification, avec ou sans planification, on 
dispose pour une sous-bande donnee W j>k de I'ensemble de ses coefficients 
d'ondelettes quantifies Q(Wj, k [n,m]). 

o On passe alors a I'etape de codage entropique dit « contextuel ». Cette 

etape se decompose en deux phases, I'une dediee a la prediction des 
esperances p et largeurs a, a partir, notamment, des premiers et seconds jeux 
de parametres de prediction determines lors de I'etape d'estimation, et I'autre 
dediee au codage entropique proprement dit. 

5 Le mode de prediction de I'esperance p et de la largeur a associees a 

un coefficient d'ondelettes quantifie depend de son contexte, ou en d'autres 
termes des coefficients d'ondelettes voisins (lorsqu'ils existent). Quatre 
situations peuvent etre rencontrees. 

En I'absence de coefficient d'ondelettes voisin nord (W[n,m-1]), c'est-a- 
- dire lorsque le coefficient d'ondelettes (W[n,1]) appartient a la premiere" ligne de 
la sous-bande concernee, on identifie les valeurs de prediction d'esperance p 
et de largeur a de la fonction laplacienne (P) aux valeurs de prediction 
d'esperance ouest /y 0 et de largeur ouest a 0 , soit p = p Q et a = a Q . 

En I'absence de coefficient d'ondelettes voisin ouest (W[n-1,m]), c'est- 
a-dire lorsque le coefficient d'ondelettes (W[1,m]) appartient a la premiere 
colonne de la sous-bande concernee, on identifie les valeurs de prediction 
d'esperance p et de largeur a de la fonction laplacienne (P) aux valeurs de 
prediction d'esperance nord p H et de largeur nord cr N , soit// =p N eta= a N . 

En I'absence de coefficients d'ondelettes voisins ouest (W[n-1,m]) et 
nord (Wtn.m-I]), c'est-a-dire lorsque le coefficient d'ondelettes (W[1,1J) 
appartient a la premiere ligne et a la premiere colonne de la sous-bande 
concernee, on effectue un codage laplacien fixe, c'est-a-dire un codage 
laplacien d'ordre 0 fonde sur la prise en compte d'une densite de probability 
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non conditionnelle. 

Enfin, en presence de coefficients d'ondelettes voisins ouest et nord, 
on identifie les valeurs de prediction d'esperance // et de largeur a de la 
fonction laplacienne (P) respectivement a la demi-somme des valeurs de 
5 prediction d'esperance nord /y N et ouest jj 0 et a la demi-somme des valeurs de 
prediction de largeur nord a H et ouest a 0 , soit/y = [(j/ N + // Q ) / 2] et a = [(a N + 
o 0 )/2\. 

Les valeurs de jj h , fj 0 , a N et a Q sont determinees a partir des relations 
(1 ) et (2) qui sont mentionnees de nouveau ci-apres pour les predictions nord 
10 et ouest : 

// N = e N .Q(W[n,m-1]) (1a) 

//o = 0o.Q(W[n-1,m]) (1b) 

o-n = an + J3 N . | Q(W[n ,m-1 ]) | (2a) 

o 0 = a 0 + 0 O . | Q(W[n-1 ,m]) | (2b) 
is ou Q(W[n,m-1]) et Q(W[n-1 ,m]) sont les coefficients d'ondelettes quantifies des 
voisins nord et ouest d'un coefficient d'ondelettes quantifie Q(W[n,m]). 

II est important de noter que le calcul de la largeur a est generalement 
biaisee lorsque la valeur du coefficient d'ondelettes quantifi6 Q(W[n,m-1]) ou 
Q(W[n-1 ,m]) associe est nulle. 

20 Dans ce cas, en constatant que Q(W[n,m]) represente I'esp6rance 

E(W[n,m] | Q(W[n,m])) de W[n,m] connaissant sa valeur quantifiee, on peut en 
deduire la relation suivante lorsque Q(W[n,m]) = 0 : 

E( | W[n,m] | | Q(W[n,m])) = 112 * \ Q(W[n,m]) | 
oO T est le pas de quantification choisi a I'etape de quantification par le 

25 quantificateur (suppose ici etre du type dit « a double boite 0 »). 

Et, partant de cette deduction, on peut alors remplacer dans les 
relations (1a), (1b), (2a) et (2b) Q(W[n,m-1]) et Q(W[n-1 ,m]), lorsqu'ils sont 
nuls, par une constante fT choisie. f est une constante strictement superieure 
a zero et inferieure ou egale a 1 . Une valeur de f egale a 1 / 2 est en theorie 

30 optimale, mais, dans le cas d'une quantification a haute resolution et de 
fonctions laplaciennes tres pentues dans la boTte 0, il est preferable de choisir f 
egale a 1/3. 
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Une fois que Ton est en possession de I'esperance p et de la largeur a 
associees a chaque coefficient d'ondelettes quantifie d'une sous-bande, on 
peut passer a la phase de codage entropique. 

Pour ce faire, il suffit de discretiser une laplacienne qui n'est pas 
s necessairement centree dans la boTte 0 du quantificateur. La discretisation 
consiste ici a echantillonner la laplacienne continue. 

La decompression de donnees d'images numeriques, compressees au 
moyen du precede de compression qui vient d'etre decrit, consiste a enchaTner 
les quatre etapes dudit precede dans I'ordre inverse. Par consequent, le 
10 precede de decompression ne sera pas decrit ci-apres. 

L'invention porte egalement sur un dispositif de compression (ou 
codeur) de donnees d'images numeriques par transformee en ondelettes. Ce 
dispositif ayant pour objet de mettre en ceuvre le precede de compression 
decrit ci-avant, seuls ses modules fonctionnels seront presentes ci-apres en 
is reference a ('unique figure. 

Le dispositif de compression D, selon l'invention, comprend un module 
de traitement PM decomposable en au moins quatre modules. 

Un premier module M1 est alimente en donnees d'image numeriques 
l[n,m] definissant des images I et est charge, d'une part, de decomposer 

20 chaque image I sur une base d'ondelettes. Cette transformation en sous- 
bandes W jik (dites LH, HL, HH et LL comme indique sur la figure 2) et sur 
plusieurs niveaux de resolution j. Chaque sous-bande est representee par un 
tableau [n,m] de coefficients d'ondelettes W jik [n,m]. II s'agit d'un module 
traditionnel, parfaitement connu de I'homme de I'art, mettant en ceuvre, par 

25 exemple, la transformee en ondelettes bi-orthogonales 7-9 de 
Daubechies/Antonini. 

Un deuxieme module M2 est alimente en coefficients d'ondelettes 
Wj, k [n,m] par le premier module de transformation Ml II est charge d'estimer, 
pour chaque coefficient d'ondelettes Wj.ktn.m] de chaque sous-bande 
30 d'ondelettes W j>k d'une image I, des premier (0 N , a N et fi H ) et second (e Q , a Q et 
fio) jeux de parametres de prediction, associes respectivement a des directions 
nord-sud et ouest-est, en fonction des valeurs des coefficients d'ondelettes de 
ses voisins nord W jik [n,m-1] ou ouest W jik [n-1,m]. Pour ce faire, le deuxieme 
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module d'estimation M2 met en ceuvre par exemple ia technique d'estimation a 
double regression decrite precedemment. 

Un troisieme module M3 est alimente en coefficients d'ondelettes 
Wj,k[n,m] par le premier module de transformation M1 et est charge de les 
5 quantifier pour delivrer des coefficients d'ondelettes quantifies Q(W j(k [n I m]). II 
s'agit d'un module de quantification scalaire (ou quantificateur) traditionnel, 
parfaitement connu de I'homme de Tart. 

Un quatrieme module M4 est alimente, d'une part, en coefficients 
d'ondelettes quantifies Q(W,, k [n,m]) par le troisieme module de quantification 

10 M3, et d'autre part, en premier et second jeux de parametres de prediction (0, 
a et /?) par le deuxieme module d'estimation M2. II est charge de coder de 
facon entropique et contextuelle les coefficients d'ondelettes quantifies 
Q(W jik [n,m]) a partir des valeurs des parametres de prediction nord (0 N , a N et 
£n) et ouest (0 O) ao et 06) des premier et second jeux. 

15 Plus precisement, le quatrieme module d'estimation M4 est agence de 

maniere a determiner pour chaque coefficient d'ondelettes quantifie 
Q(W j>k [n,m]), les valeurs de prediction d'esperance p et de largeur a de la 
fonction laplacienne P associee, compte tenu des valeurs des coefficients 
d'ondelettes quantifies de ses voisins nord Wj, k [n,m-1] et/ou ouest W jik [n-1 ,m] et 

20 des parametres de prediction nord (0 N , a N et /? N ) et ouest (0 O , a 0 et 0 O ) 
associes, puis a coder de facon entropique chaque coefficient d'ondelettes 
quantifie Q(W jik [n,m]) compte tenu des valeurs de prediction d'esperance jj et 
de largeur a qu'il a determinees. Pour ce faire, le quatrieme module de codage 
entropique M4 met en oeuvre par exemple les techniques de determination de 

25 prediction et de codage entropique contextuel decrites precedemment. 

Une fois le codage entropique termine, le quatrieme module de 
codage entropique M4 delivre sur sa sortie des donnees d'images numeriques 
compressees DC. 

Le dispositif de compression D peut egalement comporter un 
30 cinquieme module M5 charge de la planification de debit et/ou de I'allocation de 
debit par sous-bandes et/ou de la planification de qualite, preferentiellement 
selon le precede de planification presente ci-avant. Ce cinquieme module M5 
est couple au deuxieme module d'estimation M2 qui I'alimente en premier (0 N , 
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a N et £ N ) et second (6 0 , a 0 et /3 Q ) jeux de parametres de prediction et au 
troisieme module de quantification M3. 

L'invention porte egalement sur un dispositif de decompression (ou 
decodeur) mettant en oeuvre le precede de decompression precite. Ce 

s dispositif de decompression ne sera pas decrit ci-apres dans la mesure ou il 
est sensiblement constitue des memes modules fonctionnels que le dispositif 
de compression D, mais avec un fonctionnement inverse (le module de codage 
etant ici alimente en donnees d'images numeriques compressees et le module 
de transformation delivrant en sortie des donnees d'images numeriques 

10 decompressees definissant I'image initiate). 

L'invention porte en outre sur un compresseur/decompresseur, ou 
codeur/decodeur, ou encore codec, comprenant un dispositif de compression 
et un dispositif de decompression selon l'invention. 

L'invention offre un procede de codage particulierement efficace 

is assurant un gain de codage notable grace a I'exploitation de I'information 
mutuelle conjointement a la moderation laplacienne conditionnelle. 

Par ailleurs, grace aux techniques d'estimation des parametres de 
prediction et de determination des predictions, particulierement simples en 
termes de calcul, les structures des codeur et decodeur sont simplifies. 

20 En outre . du f ai* que I'on dispose d'un quantificateur different pour 

chacune des sous-bandes, il est done possible de proceder a une allocation 
locale. 

Enfin, du fait que la regie d'allocation prend explicitement en compte 
les caracteristiques du dispositif de codage entropique d'ordre 1 selon 
25 l'invention, ledit dispositif de compression est beaucoup plus efficace que ceux 
de I'art anterieur, notamment pour les applications en temps reel a haut debit 
d'entree (grandes images), en particulier dans le domaine spatial. Cette 
efficacite est encore renforcee du fait du faible nombre d'operations requises 
par le codage. 

30 L'invention ne se limite pas aux modes de realisation de procede de 

compression, procede de decompression, dispositif de compression, dispositif 
de decompression et codec decrits ci-avant, seulement a titre d'exemple, mais 
elle englobe toutes les variantes que pourra envisager I'homme de I'art dans le 
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cadre des revendications ci-apres. 
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REVENDICATIONS 



1. Procede de compression de donnees d'images numeriques par 
transformee en ondelettes, comprenant une etape de transformation des 

s donnees d'image numeriques en coefficients d'ondelettes repartis dans des 
sous-bandes, une etape de quantification desdits coefficients d'ondelettes, et 
une etape de codage entropique desdits coefficients d'ondelettes quantifies, 
caracterise en ce qu'il comprend, entre les etapes de transformation et de 
quantification, une etape d'estimation, pour chaque coefficient de chaque sous- 

10 bande de I'image, de premier et second jeux de parametres de prediction, 
associes respectivement a des directions dites « nord-sud » et « ouest-est », 
en fonction des valeurs des coefficients d'ondelettes de ses voisins nord ou 
ouest, et en ce qu'a I'etape de codage entropique, d'une part, on determine 
pour chaque coefficient de chaque sous-bande de I'image une prediction de 

is I'esperance et de la largeur d'une fonction de type laplacienne, representative 
de sa densite de probabilite, en fonction desdits premier et/ou second jeux de 
parametres de prediction et du coefficient d'ondelettes quantifie de son voisin 
nord ou de son voisin ouest, et d'autre part, on code de facon entropique 
lesdits coefficients d'ondelettes quantifies a I'aide desdites esperances et 

20 largeurs associees determines. 

2. Procede selon la revendication 1, caracterise en ce qu'a I'etape de 
codage entropique on determine pour la fonction laplacienne de chaque 
coefficient d'ondelettes quantifie, lorsque cela est possible, des predictions 
d'esperance nord et ouest et des predictions de largeur nord et ouest, ladite 

25 prediction d'esperance nord, respectivement ouest, etant definie par le produit 
entre un premier parametre de prediction du premier, respectivement second, 
jeu et la valeur dudit coefficient d'ondelettes quantifie voisin nord, 
respectivement ouest, et ladite prediction de largeur nord, respectivement 
ouest, etant definie par la somme d'un deuxieme parametre de prediction du 

so premier, respectivement second, jeu et du produit d'un troisieme parametre de 
prediction du premier, respectivement second, jeu et de la valeur absolue dudit 
coefficient d'ondelettes quantifie voisin nord, respectivement ouest. 

3. Procede selon la revendication 2, caracterise en ce qu'a I'etape de 
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codage entropique, d'une premise part, en I'absence de coefficient voisin nord 
on identifie les predictions d'esperance et de largeur de !a fonction laplacienne 
aux predictions d'esperance ouest et de largeur ouest, d'une deuxieme part, en 
I'absence de coefficient voisin ouest on identifie les predictions d'esperance et 
s de largeur de la fonction laplacienne aux predictions d'esperance nord et de 
largeur nord, d'une troisieme part, en I'absence de coefficients voisins ouest et 
nord on effectue un codage laplacien fixe, et d'une quatrieme part, en presence 
de coefficients voisins ouest et nord on identifie les predictions d'esperance et 
de largeur de la fonction laplacienne respectivement a la demi-somme des 
io predictions d'esperance nord et ouest et a la demi-somme des predictions de 
largeur nord et ouest. 

4. Precede selon I'une des revendications 1 a 3, caracterise en ce que 
I'on estime les parametres de prediction desdits premier et second jeux par 
regression. 

is 5. Precede selon la combinaison des revendications 2 et 4, caracterise 

en ce que I'on estime lesdits premiers parametres de prediction des premier et 
second jeux par une premiere regression obtenue par resolution d'une 
premiere equation W[n+1,m] = e.W[n,m) ou W[n,m] est la valeur d'un 
coefficient d'ondelettes et W[n+1,m] est la valeur du coefficient d'ondelettes 

20 voisin de W[n,m], et 6 est un premier parametre de prediction. 

6. Procede selon la revendication 5, caracterise en ce que I'on resout 
ladite premiere equation par la technique des moindres carres, ledit premier 
parametre de prediction d'un premier ou second jeu etant egal au rapport entre 
les covariance et variance non normalises des valeurs des coefficients 

2 s d'ondelettes, associes au premier ou second jeu. 

7. Procede selon la revendications 2 en combinaison avec I'une des 
revendications 4 a 6, caracterise en ce que I'on estime lesdits deuxieme et 
troisieme parametres de prediction du premier jeu par une seconde regression 
obtenue par resolution d'un premier systeme d'equations | W[n,m] - 0 N .W[n,m- 

30 1] | = a N + /? N |W[n,m-1] | , ou W[n,m] est la valeur d'un coefficient d'ondelettes, 
W ninv1 est la valeur du coefficient d'ondelettes voisin nord de W n , m , et 0 N , a H et 
/?n sont respectivement des premier, deuxieme et troisieme parametres de 
prediction du premier jeu, et en ce que I'on estime lesdits deuxieme et 
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troisieme parametres de prediction du second jeu par une seconde regression 
obtenue par resolution d'un second systeme d'equations |W[n,m] - 0 o .W[n- 
1 ,m] | ~ a 0 + fa | W[n-1 ,m] | , oCi Wn. 1>m est la valeur du coefficient d'ondelettes 
voisin ouest de W n . m , et e 0 , a Q et /3 0 sont respectivement des premier, 
s deuxieme et troisieme parametres de prediction du second jeu. 

8. Procede selon la revendication 7, caracterise en ce que I'on resout 
lesdits premier et second systemes d'equation par la technique des moindres 
carres. 

9. Procede selon I'une des revendications 1 a 8, caracterise en ce qu'en 
io presence d'un coefficient d'ondelettes quantifie voisin de valeur nulle, on le 

substitue par un produit entre le pas de quantification, utilise lors de I'etape de 
quantification, et une constante choisie, strictement superieure a zero et 
inferieure ou egale a un (1). 

10. Procede selon la revendication 9, caracterise en ce que ladite 
is constante est choisie egale a 1/3. 

11. Procede selon I'une des revendications 1 a 10, caracterise en ce qu'a 
I'etape de codage entropique, apres avoir effectue les predictions desdites 
esperance et largeur, on procede a la discretisation de la fonction iaplacienne. 

12. Procede selon Tune des revendications 1 a 10, caracterise en ce qu'a 
20 I'etape de codage entropique, apres avoir effectue les predictions desdites 

esperance et largeur, on procede audit codage par iteration. 

13. Dispositif (D) de compression de donnees d'images numeriques par 
transformee en ondelettes, comprenant des moyens de traitement (PM) munis, 
d'une premiere part, d'un module (M1) de transformation de donnees d'image 

25 numeriques en coefficients d'ondelettes repartis dans des sous-bandes, d'une 
deuxieme part, d'un module (M3) de quantification des coefficients 
d'ondelettes, et d'une troisieme part, d'un module (M4) de codage entropique 
des coefficients d'ondelettes quantifies, caracterise en ce que lesdits moyens 
de traitement (PM) comprennent en outre un module d'estimation (M2) agence 

so pour estimer pour chaque coefficient de chaque sous-bande de I'image, de 
premier et second jeux de parametres de prediction, associes respectivement 
a des directions dites « nord-sud » et « ouest-est », en fonction des valeurs des 
coefficients d'ondelettes de ses voisins nord ou ouest, et en ce que ledit 
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module de codage entropique (M4) est agence, d'une part, pour determiner 
pour chaque coefficient de chaque sous-bande de I'image des valeurs de 
prediction de I'esperance et de la largeur d'une fonction de type laplacienne, 
representative de sa densite de probability en fonction desdits premier et/ou 
s second jeux de parametres de prediction et du coefficient d'ondelettes quantifies 
de son voisin nord ou de son voisin ouest, et d'autre part, de coder de fagon 
entropique lesdits coefficients d'ondelettes quantifies a I'aide desdites 
esperances et largeurs associees determinees. 

14. Dispositif selon la revendication 13, caracterise en ce que ledit 
lo module de codage entropique (M4) est agence pour determiner pour la 

fonction laplacienne de chaque coefficient d'ondelettes quantifie, lorsque cela 
est possible, des valeurs de prediction d'esperance nord et ouest et des 
valeurs de prediction de largeur nord et ouest, ladite valeur de prediction 
d'esperance nord, respectivement ouest, etant definie par le produit entre un 

is premier parametre de prediction du premier, respectivement second, jeu et la 
valeur dudit coefficient d'ondelettes quantifie voisin nord, respectivement ouest, 
et ladite valeur de prediction de largeur nord, respectivement ouest, etant 
definie par la somme d'un deuxieme parametre de prediction du premier, 
respectivement second, jeu et du produit d'un troisieme parametre de 

2 o prediction du premier, respectivement second, jeu et de la valeur absolue dudit 
coefficient d'ondelettes quantifie voisin nord, respectivement ouest. 

15. Dispositif selon la revendication 14, caracterise en ce que ledit 
module de codage entropique (M4) est agence, d'une premiere part, en 
I'absence de coefficient voisin nord pour identifier les valeurs de prediction 

is d'esperance et de largeur de la fonction laplacienne aux valeurs de prediction 
d'esperance ouest et de largeur ouest, d'une deuxieme part, en I'absence de 
coefficient voisin ouest pour identifier les valeurs de prediction d'esperance et 
de largeur de la fonction laplacienne aux valeurs de prediction d'esperance 
nord et de largeur nord, d'une troisieme part, en I'absence de coefficients 

jo voisins ouest et nord pour effectuer un codage laplacien fixe, et d'une 
quatrieme part, en presence de coefficients voisins ouest et nord pour identifier 
les valeurs de prediction d'esperance et de largeur de la fonction laplacienne 
respectivement a la demi-somme des valeurs de prediction d'esperance nord 
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et ouest et a la demi-somme des valeurs de prediction de largeur nord et ouest. 

16. Dispositif selon I'une des revendications 13 a 15, caracterise en ce 
que ledit module d'estimation (M2) est agence pour estimer les parametres de 
prediction desdits premier et second jeux par regression, 
s 17. Dispositif selon la combinaison des revendications 14 et 16, 

caracterise en ce que ledit module d'estimation (M2) est agence pour estimer 
lesdits premiers parametres de prediction des premier et second jeux par une 
premiere regression obtenue par resolution d'une premiere equation W[n+1,m] 
= 0.W[n,m] oCi W[n,m] est la valeur d'un coefficient d'ondelettes et W[n+1,m] 
io est la valeur du coefficient d'ondelettes voisin de W[n,m], et 6 est un premier 
parametre de prediction. 

18. Dispositif selon la revendication 17, caracterise en ce que ledit 
module d'estimation (M2) est agence pour resoudre ladite premiere equation 
par la technique des moindres carres, ledit premier parametre de prediction 

is d'un premier ou second jeu etant egal au rapport entre les covariance et 
variance non normalises des valeurs des coefficients d'ondelettes, associes 
au premier ou second jeu. 

19. Dispositif selon la revendications 14 en combinaison avec I'une des 
revendications 16 a 18, caracterise en ce que ledit module d'estimation (M2) 

20 est agence, d'une part, pour estimer lesdits deuxieme et troisieme parametres 
de prediction du premier jeu par une seconde regression obtenue par 
resolution d'un premier systeme d'equations |W[n,m] - 0 N .W[n,m-1]| = a N + 
0h | W[n,m-1] | , ou W[n,m] est la valeur d'un coefficient d'ondelettes, est 
la valeur du coefficient d'ondelettes voisin nord de W n>m , et 0 N> a N et fi H sont 

25 respectivement des premier, deuxieme et troisieme parametres de prediction 
du premier jeu, et d'autre part, pour estimer lesdits deuxieme et troisieme 
parametres de prediction du second jeu par une seconde regression obtenue 
par resolution d'un second systeme d'equations |W[n,m] - 0 o .W[n-1,mj | = a Q 
+ 0o | W[n-1,m] | , ou W n . 1>m est la valeur du coefficient d'ondelettes voisin ouest 

so de W„, mi et 9 0 , a 0 et p 0 sont respectivement des premier, deuxieme et 
troisieme parametres de prediction du second jeu. 

20. Dispositif selon la revendication 19, caracterise en ce que ledit 
module d'estimation (M2) est agence pour resoudre lesdits premier et second 
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systemes d'equation par la technique des moindres carres. 

21. Dispositif selon Tune des revendications 13 a 20, caracterise en ce 
que ledit module de codage entropique (M4) est agence, en presence d'un 
coefficient d'ondelettes quantifie voisin de valeur nulle, pour le substituer par le 

s resultat d'un produit entre le pas de quantification, utilise par le module de 
quantification, et une constante choisie, strictement superieure a zero et 
inferieure ou egale a 1 . 

22. Dispositif selon la revendication 21, caracterise en ce que ladite 
constante est choisie egale a 1/3. 

io 23. Dispositif selon Tune des revendications 13 a 22, caracterise en ce 

que ledit module de codage entropique (M4) est agence pour proceder a une 
discretisation de la fonction laplacienne apres avoir effectue les predictions 
desdites esperance et largeur. 

24. Dispositif selon Tune des revendications 13 a 22, caracterise en ce 

15 que ledit module de codage entropique (M4) est agence pour proceder audit 

codage par iteration apres avoir effectue les predictions desdites esperance et 
largeur. 

25. Codec de donnees d'image numeriques, caracterise en ce qu'il 
comprend un dispositif de compression (D) selon I'une des revendications 13 a 

20 24. 

26. Procede selon la revendication 1, caracterise en ce qu'il permet une 
regulation de debit en temps reel par la mise en ceuvre d'un algorithme de 
planification de debit dit "avec allocation de debit" consistant a laisser un 
module de planification, d'une part, choisir independamment les pas de 

25 quantification, utilises lors de la quantification, des differentes sous-bandes de 
sorte que le debit total soit egal a une consigne de debit, et d'autre part, choisir 
la facon dont ce debit doit etre reparti dans les differentes sous-bandes. 

27. Procede selon la revendication 1, caracterise en ce qu'il permet une 
regulation de debit en temps reel par la mise en oeuvre d'un algorithme de 
planification de debit dit "avec allocation de debit pure" consistant a lier les pas de 
quantification des differentes sous-bandes, en general identiques, un module de 
planification se contentant de choisir le pas de quantification unique pour toutes les 
sous-bandes qui permettra d'atteindre une consigne de debit fixee. 
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